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Abstract 

We explore some aspects of "reconstructing" the heavy singlet sector of supersymmetric type I seesaw models, 
for two, three or four singlets. We work in the limit where one light neutrino is massless. In an ideal world, 
where selected coefficients of the TeV-scale effective Lagrangian could be measured with arbitrary accuracy, the 
two-singlet case can be reconstructed, two three or more singlets can be differentiated, and an inverse seesaw with 
four singlets can be reconstructed. In a more realistic world, we estimate la — >■ ipy expectations with a "Minimal- 
Flavour- Violation-like" ansatz, which gives a relation between ratios of the three branching ratios. The two singlet 
. model predicts a discrete set of ratios. 
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^ 1 Introduction 

New particles which are too heavy to be produced on-shell, can nonetheless leave observable traces of their virtual 
exchange. At scales well below their mass, their effects can be described by an effective Lagrangian containing non- 
CIh" renormalisable operators induced by the exchange of the heavy new particles. While it is clear that the effective 
^ Lagrangian can always be constructed from a New Physics model, the prospects for "reconstructing" the New Physics 
from the effective Lagrangian are nebulous. Such a "reconstruction" would be interesting for any New Physics at a 
. scale beyond the reach of the LHC. 

This paper aims to analytically explore this reconstruction in the simplest of toy models. We suppose that light 
. neutrino masses, with one massless neutrino, are generated by a supersymmetric type I seesaw model[l], with two, 
' three or four heavy singlets of mass Mj ^ TeV ^ . When the number of singlets njv does not exceed the number of 
lepton doublets, this model is in principle "reconstructable" [3], that is, the masses and mixing angles of the singlet 
'■ sector can be computed from parameters (masses and mixing angles) of the doublet sector. This reconstruction has 
QQ , been studied from numerous perspectives [4, 5, 6, 7, 8, 9]. Analytic formulae for the reconstruction of two-singlet 
models were presented in [5]. We wish to know what we could learn about the singlets, in principle and also with 
: some degree of realism, from effective operators involving doublets. So we start in section 3 with a model containing 
" ^ I ' two singlets, which has sufficiently few parameters that it could be disfavoured by observations. Then in section 4 
^ . we consider models with three singlets, in two limits: where there is a massless doublet neutrino because a Yukawas 
eigenvalue vanishes, and the case where a singlet is infinitely massive. As expected, the three-singlet sector is in 
principle reconstructable when there are three non-zero Yukawa eigenvalues, but not when there are only two, because 
in this latter case the singlet components which couple via the zero Yukawa are decoupled. Finally, in section 5, by 
considering even more exotic operators, we show that in principle it could be possible to distinguish between two, 
three or more singlets models, and to reconstruct a four singlet "inverse seesaw" [10] model. 
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2 Notation and Assumptions 

The seesaw model [1] is a natural and minimal extension of the Standard Model which fits the observed neutrino 
masses. We consider a supersymmetric seesaw model (with conservation of R-parity^) because we are interested in 
seesaw reconstruction; the slepton masses may contain additional information about the seesaw parameters [3]. Also, 
in the non-supersymmetric seesaw, there are 0{{y'^ AI)'^) contributions to the Higgs mass which must be fine-tuned. 
At scales above the singlet masses Mj, the superpotential can be written 

Wiep = [Y-]o.c.{Lc.Hd)E-^ + [Yn^,/(L^i/„)Aj + ^N}N}. (1) 

*E-mail address: s.davidson@ipnl.in2p3.fr 
tE-mail address: m.elmer@ipnl.in2p3.fr 

^This means we do not consider "low-scale" seesaws, with singlet masses in the eV — >■ TeV range [2]. 

^We assume that neutrino masses arise due to the seesaw, so R-parity, or some other symmetry, must prevent [11] other dangerous 
lepton number violating interactions which could generate majorana neutrino masses [12]. 
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This expression is in the eigenbases of VY^^ for the {La}, V^'V for the {E^^}, and M for the {iV/}. These choices 
correspond to the charged lepton mass basis, refered to as the flavour basis and labeUed by greek letters, and the mass 
eigenstate basis of the heavy singlets, labelled by roman capitals which run from 1 to Un- The doublet contraction 
is antisymmetric (LHd) = ElH° — NlH~ , the Yukawa indices are ordered left-right, and we will allow un = 2,3,4 
generations of singlets N, whose masses M/ will usually be taken ^ TeV. The resulting Lagrangian is 



£ = y'J^'me'^ + e''H:[Y'^*UNj + ^N-iNi + ... + h.c. 



(2) 



where G {ue, Vf^jUr} are the charged lepton Yukawa couplings, the singlet neutrinos are written as four-component 
fermions, and the ... includes sparticle interactions. 

The number of parameters [13] in the superpotcntial (1), or Lagrangian (2), will depend on the number of singlets 
Uat, and on whether the elements of the matrices M and Y'' are allowed to be complex or restricted to be real. In 
the case of complex matrices with n^v singlets, there are un masses {Mj}, which can be taken real by a phase choice 
on the {Nj}, three eigenvalues {y^}, which can be taken real by a relative phase choice between the {E^} and {La}, 
and 3 x hn complex entries in Y'^, from which 3 phases can be removed by suitably choosing the phase differences 
between the three doublets {La} and the singlets {Nj}. So we expect 7 x real parameters. If the matrices M and 
Y" are restricted to be real, there would be 3 + 4 x tiat real parameters. 

The neutrino Yukawa matrix Y" is a 3 x njv matrix, with at most mm{3,njv} non-zero eigenvalues. It can be 
diagonalised by independent unitary transformations V/, (which acts in on the left and is 3 x 3) and V/j (which is 
njv X TiN and acts on the right). So in the case of njv = 3 : 



VlY^'V 



Dy- = diag {2/1,^2,2/3} 



(3) 



For njv = 2, the same formula can be used, but with yi 
corner of the 3x3 matrix appearing in eqn (3): 



0, and by putting the non- trivial 2 x 2 Vr in the lower right 



[V. 



R\3x3 = 



1 

[Vfi]2x2 





(4) 



The case where uat = 4, can be similiarly delt with, by embedding the 3 x 3 Vl into a 4 x 4 matrix. 

At scales <C M, where the N are not present as on-shell particles, the effective Lagrangian will contain non- 
renormalisable operators induced by A'^ exchange. For simplicity, in this paper we focus on CP conserving observables 
(for instance, we neglect electric dipole moments 

At dimension flve, arises the majorana mass operator for the doublet neutrinos 



(5) 



where Vu = {Hy) = vsml3 ~ w = 174 GeV, Dm = diag{mi,m2,ms}, and the neutrino mass differences are taken to 
be 



lAm^il ~ 2.4 X 10"-'' cV^ \i 



lAm^il ~ 7.6 X 10"^ eV^ (6) 

We assume that the lightest doublet neutrino is exactly massless (relaxing this assumption is discussed at the end 
of section 4.1), and consider separately the hierarchical(NH: ms > m2 > mi = 0)) and inverse hierarchical (IH: 
1712 > mi > tn^ = 0) patterns. U is the leptonic mixing matrix 



U = 



p C13C12 



S12C13 



-iS 



e'" (-S12C23 - S23S13C12 e**^) e^P (C23C12 - S23S13S12 ) 
e*" (S23S12 - S13C23C12 e'*) e'"^ (-S23C12 - S13S12C23 e"' 



Si3 e 
S23C13 
C23C13 



(7) 



where Cij = cos 0ij, etc., ^23 



7r/4[15], sin^0i2 = 0.3 [15] (or sometimes sin^ ^12 — 1/3), and Ues = 0.18, the central 
value of the recent T2K indication[16] of non-zero ^13 (averaged over NH and IH). Since a neutrino is massless, one 
of the Majorana phases a, (3 vanishes. In practise, we will usually neglect all phases for simplicity (we do not discuss 
CP violating observables). 

At dimension six, the seesaw model induces the lepton number conserving operator [17, 18] {£aH*) f) (i^H) with 
coefflcient 



[d]a^ = 



Y''M-^M-i*Y'^^ 



al3 



(8) 



After electroweak symmetry breaking, this operator modifies the neutrino kinetic terms, so can contribute to the 
non-unitarity of the leptonic mixing matrix[19]. 

The third operator which we wish to use for seesaw reconstruction is the dimension six electromagnetic dipole 
operator, which can be written, after electroweak symmetry breaking, as 



(9) 



'For a review and references, see e.g. [14] 



BR 


current bound 


future 


T — )• /i7 

r — )• 67 


2.4 X 10-i2[22] 

2.5 X lO-'^pS] 
1.9 X 10-^[23] 


~ 10"", (MEG [22]) 
~ 10"^, (super-B factories [24]) 
~ 10"^, (super-B factories [24]) 



Table 1: Current bounds and hoped-for sensitivities to lepton flavour violating branching ratios, normalised to leptonic 
weak decays, as in eqn (12). 



We assume that the coefficients [XiJ^^ of this chirality-flipping operator have the "Minimal Flavour Violation [20] -like" 
form 

Xl = vCfY''^ + ?;C^Y^^Y''Y''1' , (10) 

where v ~ 174 GcV is the Higgs vcv, and the coefficients have mass dimension -2 and are assumed flavour in- 
dependent. Information about the flavour off-diagonal elements of [VV^] therefore can be obtained. This form is 
"MFV-like", because the lowest order dimcnsionlcss lepton flavour-changing operator in the seesaw model is [Y'^Y'^^]. 
This parametrisation of is motivated by supersymmctry, where loop corrections to the soft masses can approxi- 
mately give such a dependance on [Y'''Y'''^]. In the Appendix, we review the supersymmetric motivation for eqn (10), 
and estimates for the coefficients and . These imply that 

[yV^Xe ^ 2x10-3 , (11) 

to suppress the rate for n ej. We impose this "approximate zero" on our models, because it approximately fixes 
a parameter, with a weak dependence on the vahic of . We try to avoid assuming a value for , because a 
precise calculation depends on several model parameters at various scales. This means we extract parameters from 
— or predict — ratios of lepton flavour violating rates (such as eqn (13)), but we cannot predict the rates (such as eqn 
(12)). 

Lepton flavour violating radiative decays, £a — >■ ^/ST, proceed via the operator (9) at a rate given by 

m.^h ) = »<(\^^^^\' + \^n.,f)^ . (12) 

Table 1 lists the current bounds, and hoped for sensitivities of running or planned experiments. If two different lepton 
flavour violating decay rates were observed, for instance n ^ ej and r — )• ^j, then the approximation (10) implies the 
following equality 

g(M ^ e7) ^ |[Y'-Y'-t]^^|2 
Bi?(r^^7) l[Y-Y'^^].^P 

We assume that the eigenvalues of Y-^Y-^t ^ vIdIVl are sufficiently hierarchical, that in — >■ £/3"/ processes, only 
the terms involving the largest eigenvalue ys need to be considered: 

|[Y''Y^t]^^|2_^j^4|^^^^^^*^^|2 _ (14) 

As all known Yukawa eigenvalues of the SM are hierarchical, this assumption is not unreasonable for the neutrino 
Yukawa matrix. The approximation (14) implies a relation among the three rates ^01, because they are all 

controlled by one row of Vl, which depends on two angles. The y| in eqn (14) cancels in ratios of branching ratios 
such as eqn(13). The — > £/3j branching ratios are oc y| (an allowed range for ys might be disentangled from slepton 
masses, if the spectrum of supersymmetric particles was well-known); we set j/f = 1 to estimate branching ratios. If 
the assumption of sufficiently hierarchical Yukawa eigenvalues is dropped, the £a — ^/37 rates depend on additional 
parameters. This destroys the relations between different rates. 

The aim of this study, is to "reconstruct" the singlet sector of seesaw models. With this aim, we consider three 
matrices in doublet flavour space: the light neutrino mass matrix [m^] of eqn (5), the dimension six [18] operator [d] 
of eqn (8), and the matrix 

P = = VIt>IVl (15) 

The flrst two arc coefficients of operators in the effective lagrangian, and the off-diagonal elements of P could be 
extracted from the slepton mass matrix. To finish determining Y'^Y''^ (in the approximation of neglecting phases), 
three more inputs would be required, which we take to be the eigenvalues. We assume they are hierarchical, that 
ys = 1, and wc will keep the remaining two as free parameters. We do not study the extraction of these eigenvalues 
from slepton mass matrices, although this can be envisaged [6, 5] in some circumstances. 

If the elements of the light neutrino mass matrix, and of either of P or [d] could be measured exactly, then the 
parameters of the high scale Lagrangian of eqn (2) can be reconstruct in the njv < 3 seesaw. This is clear in the case 



of P and [mi^]: diagonalising P — Y'^Y"^ gives the Y" eigenvalues and the unitary matrix Vl which transforms from 
the doublet basis where is diagonal, to the basis where Y"^ is diagonal. The matrix 

Wl=VlU , (16) 

which tranforms from the neutrino mass basis to the basis Y'^ diagonal, is therefore known. So the expression for the 
mass matrix of massive light neutrinos 

D„ = Wl^-D,VrT>m'Wr^T>,WIvI (17) 

can be inverted for Vr and = diag{l/Mi, ...}: 

VrDm-'Vr^ = T>,-^WLTimWlT>Z^vl (18) 

provided that has no zero eigenvalues. The matrix Wl parametrises the deviation of Vl (whose angles control 
LFV processes) from the leptonic mixing matrix U . If there are only two non-zero eigenvalues in Y"^, so the lightest 
neutrino is massless, then Wl is a 2 x 2 matrix, and Vl differs from U hy a, single rotation. 

As discussed by Broncano,Gavela and Jenkins[18], the parameters of the high-scale singlet lagrangian can also be 
reconstructed from [m,y] and [d]. Although [d] is the coefficient of a dimension six operator, (so suppressed oc 1/M^), 
if it was known exactly, then 

v^[ui,]-^[d] [m.*]-i = [Y-^]-i[Y'']*-i (19) 

so the matrices [Y^^] and [M] can be obtained. 

In an ideal world, where two of the the three matrices P, [m^] and [d] could be measured, then the third could be 

used as a test of the seesaw. 

In summary, we assume that the lightest doublet neutrino is massless, that ^13 ~ .18, and we neglect the phases 
of the lepton mixing matrix. We assume that the neutrino Yukawa coupling matrix Y"^ has hierarchical eigenvaluess 
with 2/3 = 1, and that flavour off-diagonal elements of [Y'^Y''^] 

can be extracted from the rates for £q y i-B^- 



3 Two singlets 

The first model that we consider contains two heavy singlet ("right-handed") majorana neutrinos Nj. This is the 
simplest see-saw model that can explain all observed neutrino oscillations, and has been carefully studied by several 
people [25, 5]. Here, we review the expectations for 1^ f-pj decays, and the prospects for reconstructing the singlet 
sector. This model predicts that one of the active doublet neutrinos has a vanishing mass; in the case of the normal 
hierarchy, this massless neutrino is taken to be ui, whereas it is U3 for the inverse hierarchy. 

In general there are two possibilities to write down this model. The first possibility is to use a 3 x 2 neutrino 
Yukawa coupling matrix [Y'^j^/, so that it can connect the 3 generation doublet lepton sector to the 2 generation 
singlet sector. The second possibility is to write also the right handed sector as 3 x 3 matrices and fill up the unused 
components with zeros, so that all matrices that are used can by diagonalised. We will use the latter one. 



3.1 Real case 

In this first step, all phases are set to zero. 

To reconstruct the high scale parameters, using equation (18), we need the eigenvalues of Y"^, and the matrix Wl, 
which parametrises the deviation of Vl from U. In the two singlet model, equation (18) can be written in terms of 
2x2 matrices, because the basis of the charged Icptons is not used. We take 1/3 = 1,2/2 remains as a free parameter, 
and writing 

[VlU = [WLhK^ (20) 

allows to obtain the orthogonal 2x2 matrix Wl from lepton flavour violating decay ratios, as given in eqn (14). The 
index j runs from 2 to 3 for the normal hierarchy, and 1 to 2 for the inverse hierarchy. 

The bound (11) from the non-observation of n ej [22], combined with the assumed hierarchy 1 = 2/3^ 2/2 (see 
eqn (14)), approximately fixes the angle 9wl, whose value hardly changes as |VL3eV^3j(i| varies from 2 x 10"'^ to zero. 
Of course this fine-tuning of an element of Vl changes the branching ratios which depend on it, but has little effect 
on the reconstructed inverse mass matrix of the singlet neutrinos. The "approximate zero" of eqn (11) can be 

implemented either as |VL3e| = \swLUej + cwlUsj+iI < 10"^ or |Vl3/x| = IswlU^^j + cwLU^j+i \ < 10"^ (where j = 1 
for the inverse hierarchy, j = 2 for a normal hierarchy, and sw = smOwL, etc.). These two possibilities determine the 
two possible solutions for Owl- If \VL3e\ is small, then the branching ratio for r — )• /i7 could be detectable, and 

BRjl^^e^) ^ \Vl3,\' ^ \~U:j+iU^j+U:jU,j+i? 
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Whereas if |VL3/i| is smah, then the branching ratio for t — > 67 could be detectable, and 

BRifl ^ 67) \VL3e\' I - C/;, + lC/e, + U^^U^j + l ^ 



br{t ^ fij) \VL3r\' I - u;^^,Urj + u;^Urj+i\ 



0{1) (22) 



We now consider the expectations for these ratios of branching ratios, for the normal and inverse hierarchy patterns 
of light masses. 

For the inverse hierarchy, the case Vlsb gives 

S^i^ = 1 (23) 
BR{t 67) 

which implies that r — > 57 is beyond the reach of planned experiments, see table 1. However, t — > jij oc jVs^Vg*^!^ ~ .25, 
so according to the SUSY estimate for (see eqn (56)), BR{t — > /X7) ~ 2.3 x 10"'', so should be seen in upcoming 
experiments (see table 1). 

The case Vts^ ^ in the inverse hierarchy is curious, because Vl3i- — Si3 and Vlzc 1- Therefore r ^ /i7 is 
beyond the reach of planned experiments : 

- si, , (24) 

BRin^ej) 

and since IVtsrV^sel — sis, the supersymmetric estimate of eqn (56) gives BR{t 67) ~ 10~^s1,. This could be 
detectable, for cooperative supersymmetric masses and mixing angles, and S13 close to the current T2K central value 
~ .18. Unfortunately, with uncooperative supersymmetric parameters, all — > ^,97 decays could be out of reach. 

Consider now the case of a normal hierarchy for the light neutrino masses, so the Vtsa are defined with j = 2.. 3 
in eqn (20). If |VL3e| — 0, then we obtain 

g(/Z ^ 67) ^ I - Ue3U^2 + Ue2U^3\' ^ ^ 3^ 

which implies that t — > 67 is beyond the reach of planned experiments. However, the supersymmetric estimate for 
from the Appendix implies that BR{t — > /U7) ~ 10"''. Alternatively, when |Vt3^| — >■ 0, one obtains BR{t — >■ 67) ~ 
1 X 10~^, and BR{t — >■ /X7) out of reach at a branching ratio of order BR{^ — >• 67): 

BRjr ^ M7) ^ I - UUUr2 + t/^^t/.aP ^ 4 

Notice that in the normal hierarchy case, it is not possible to suppress all the £a ^/37 decays, as was possible for 
the inverse hierarchy: there is always one small branching ratio of the order of fj, ^ 67, and one BR{t — )• t'j) ~ 10~^, 
which could be measurable at a superB factory. 

We can now attempt to reconstruct the mass of the heavy singlet neutrinos, with 9wl fixed by eqn (11). The 
masses are similar for the four possible values of 9wl- The mass M2 of the lighter singlet depends on where y2 is 
the yet unfixed second neutrino Yukawa eigenvalue. Thermal leptogenesis [26] ^ prefers M2 ^ 10^° GeV, which can 
be obtained for 

0,002 < j/2 < 0,007 normal hierarchy (27) 
0, 003 < y2 < 0, 009 inverse hierarchy. (28) 

In both cases, the decay rate r2 of this lighter singlet is ~ 45i?, where H is the Hubble expansion rate at T M2, 
so leptogenesis would take place in the strong washout regime. In these limits the heavier singlet mass is found to be 
M3 « 2 • 10^^ GeV in the hierarchical case and M3 « 6 • 10^^ GeV in the inverse hierarchy. 

3.2 complex case 

In this section, we check that the (unknown) phases of Wl can always be chosen to compensate those of U, so as to 

obtain negligeable imaginary parts on the matrix Vl- We perform this check for a hierarchical light neutrino spectrum. 
Recall that the matrix elements of Vl control the £a — ^ ^^7 rates. If we could not arrange their imaginary parts to be 
negligeable, then our approximation of neglecting phases would not work. 

As only two of the light neutrinos have masses, only one Majorana phase is required in the lepton mixing matrix 
U, so we set a = 0. The 2x2 unitary matrix Wl can be written 

/I 

W = 6^"^" I cwlb'^^^ -swhe-'"!"^ \ . (29) 




*Our model is supersymmetric, so some solution [27] to the gravitino problem would be required. 
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As in the previous section, wc assume that 



J^g|i(7+a+2;3-25)^^^^^^ + Ci3S12SIVl) • {cwLCliS2S + ei'^^+^+^^hwL{ci2C23 - e-'^Si2Sl3S23)) | (30) 



must be small to suppress /x — >■ 67, from which we want to determine the mixing angle 9wl- To obtain a value of the 
order of 10""^ or even smaller, both the real and imaginary parts of one of the two parentheses must be small. The 
first parenthese is small for 



^(7 + (7 + 2/3 - 2(5) « 0, ±7r, ±2tt, ... 

Owl ~ arctan ( ) = 0.316 

VC13S12/ 



(31) 



and the second is small for 



1/ noN /cos (5- /cos5-11.2 

-(7 + CT + 2/3) w arccot . ^^^^^^^^^ = arccot 



2 ^ V sin ^ J \ sin 5 y (32) 

ewL « -0.77 



We see that the condition of small /i 67 branching ratio restricts the possible phases. 

An important question is whether the predictions of the other two branching ratios depend strongly on the choice 
of the phases. We find that these two decay rates also only depend on the sum of 7 + ct. If the phases vary in such a 
way that they still satisfy equation (31) or (32) and if |Vl3^jF£3,, | stays constant, we find that the other two branching 
ratios vary. But their overall variations over the whole parameter space are not bigger than a factor two. This factor 
two might be crucial for an actual determination of the branching ratios, but in the scope of our approximations this 
factor two will not influence our conchisions. 

The phases also have an influence on the reconstructed Majorana masses. The first observation that we make is, 
that the Majorana mass eigenvalues become complex. However by one phase rotation we can absorb one phase in 
the W matrix. So finally we only have one complex mass eigenvalue. Finally we look at the absolute values of our 
Majorana masses and compare them to their value in the real case. We find that the variation does not exceed a factor 
four. This will also not have any influence on our conclusions. 

In the following we will neglect phases in our reconstruction procedure. As we are not interested in CP violation 
and the variations of the physical observables for different phases are rather small, we will place all our following 
analysis in a real case. 



4 Three singlets 

This section studies the prospects for reconstruction of a seesaw model with three heavy singlets Nj. In this model, 
all the light neutrinos can be massive. We study the limit where the smallest light neutrino mass m,^,min 0, to 
represent the experimental situation where nij^^min \/ ^"^soi unknown. This limit could arise [5] when one of 
the singlets is very heavy, I/M3 0, so does not contribute to the observed light masses, or when one of the 'Y" 
eigenvalues docs not contribute to the light masses, y\ 0. We wish to know if such three singlet models can be 
distinguished from the two-singlet model of the previous section. All phases are neglected (set to zero) in this section. 

4.1 One infinite mass eigenvalue 

As regards the light neutrino mass matrix [m], the three singlet case with I/M3 — > is similar to the model of two 
singlets, because makes no contribution to [m] by construction. If, in addition, makes no contribution to the 
coefficients Xp of the electromagnetic dipole operator (see eqn (33)), then iVs is decoupled and the model looks like 
a two-singlet model. 

The aim here is to consider the situation where N-^ did contribute to Xp, but not to [m]. In a supersymmetric 
scenario, this situation could arise in the very narrow window where M3 < Mx (recall that Mx ^ mpi is the scale 
at which the soft masses appear), but M3 is irrelevant for [m]. It is clear that when the mass M3 exceeds the scale 
Mx, then the heaviest singlet is not present to contribute to loop corrections to the soft masses. The combination of 
Yukawas in eqn (10), would therefore be replaced by 

where P12 is the projector onto the two lighter singlets, in the Nj mass eigenstate basis. It is straightforward to check 
that this matrix has two non-zero eigenvalues. If it is identified as a 2 x 2 matrix Y'^Y'^^ , then, in combination with 



vlB,VRPi2V^f,B^VL 



12 



1 
1 




(33) 



the light neutrino mass matrix [m], it allows to reconstruct the seesaw model of the two lighter singlets. The heaviest 
is fully decoupled and irrelevant. 

For the remainder of this subsection, we assume that the Yukawa matrix Y'' has three non zero eigenvalues, so 
that Wl (or Vl) and Vu are 3x3 matrices, described by three mixing angles. If Vl is written in the form of the 
leptonic mixing matrix U, then with the approximation (14), two of the angles could be determined from two ratios 
of branching ratios : 

BRjf^ ^ ^ I^L3j^ ^ 4^23 BR{p ^ e7) ^ |l^L3eP ^ 4l13 (34) 

That means that measuring one of these ratios of branching ratios, does not predict the other (This is different from the 
two singlet model of the previous section). However, if | Vts/^V^ggl is fixed by an observation of /x — )■ e-y, then measuring 
one ratio of eqn (34) determines the other — this is illustrated in figure 1, where the upper diagonal (black) corresponds 
to BR{fi — > 67) = 2.4 X 10~^^, and the lower diagonal (red) line would correspond to BR{^ — >■ 67) = 3.0 x 10~^^ 
(The branching ratios are estimated via eqn (56)). This correlation is not a consequence of the supersymmetric seesaw 
model; rather, it follows from assuming that: 

1. the coefficient-matrix of the dipole operator contains one lepton flavour changing matrix as in eqn (10), 
which is dominated by its largest eigenvalue, and, 

2. that the small branching ratio of /i — > 67 imposes an approximate zero in X^. 

So if the prediction of figure 1 was verified, it would be consistent with the three-singlet seesaw, as well as any other 
MFV-like model satisfying the above two assumptions. And if some other ratios of branching ratios were measured, 
it could likely be accomodated in a supersymmetric three-singlet seesaw model by relaxing the above assumptions. 




1E-14 1E-13 1E-12 1E-11 1E-10 1E-9 1E-8 1E-7 



br{t - ^n ) 



Figure 1: Branching ratios for t — > /i7 and r ^ 67 (normalised to leptonic decays as in eqn (12)) in the hypothetical 
cases that fi —>■ was observed, with BR{fi — >■ 67) = 2.4 x 10~^^ (upper rising-towards-the-left black line), or 
BR{^ ej) — 3.0 X 10"^'^ (lower rising-towards-the-left red line). The correlation arises from the assumptions 
discussed after eqn (34). The BRs are estimated with eqn (56). The estimated sensitivity of superB factories is 
outside the (blue) rectangle. This plot can also be interpreted as an upper bound on the BRs which are consistent 
with a. fj, ej bound: the allowed region is to the left and below the diagonal line sloping- up-to- the- left. The 
horizontal and vertical cutoffs are uncertain because they are determined by our estimate of . 

To reconstruct the singlet masses with equation(18), requires the neutrino mass matrix uii, and the whole P matrix. 
However, in the approximation (14), neither the Y" eigenvalues nor the third mixing angle of Vl can be determined 
from rare decays (degenerate sneutrinos [28, 6] could give additional information). Our reconstructed singlet masses 
therefore depend on the unknown 0yLi2, yi and 2/2 ■ We can estimate ranges for these parameters which give Mi ~ 10^'^ 
GeV, as would be suitable for thermal leptogenesis. 

The inverse mass matrix of the singlets is given in equation(18), as a function of the "known" doublet parameters. 
When one eigenvalue of a 3 x 3 matrix M vanishes, the two remaining eigenvalues are given by [29] 



Mi,M2 = i 



Tr[M] ± ./(Tr[M])2 - 4(AfiiAf22 + M11M33 + Af22M33 - A^ls - ^ M^z) 



(35) 
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We find that for generic choices of the mixing angles in Vl, the singlet masses Mi and M2 have the form ~ v^yl/m, 
as plotted in figure 2, where m — -^/Am^ + -y/Aml^. This is unsurprising, since the light masses corresponding to 
eqn (6) are much more degenerate than the squared eigenvalues of . So a given M/ is usually controlled by a single 
yf, and the singlet mixing angles in Vr are generically small. That is, Vr usually can be written as a rotation matrix 
with small mixing angles multiplied by either the identity matrix or a permutation. We find that we can always choose 
Yukawa eigenvalues to obtain Mi ~ 10^° GeV, as suggested for thermal leptogenesis. This tells us that even if all rare 
decays were measured, the lightest singlet mass is not determined, so whether it has a value suitable for leptogenesis 
is not known. 

To reproduce, in this model, the two singlet case of section 3, can be achieved by setting 9wli2 and OwLis in Wl 
to zero. Then yi becomes irrelevant for the reconstruction, 0vfl23 plays the role of Owl in the two singlet model. In 
this limit. Mi depends quadratic on y2. 




0,01 0,1 



yi 



Figure 2: General dependence, in the three-singlet model where I/A/3 — > 0, of the singlet masses on the smallest 
Yukawa eigenvalue yi. In this plot, j/2 — -01 and j/3 — 1. The angles in Wl were taken similar to those of U: 
OwLi3 = -18, 0wLi2 — -6, and 9wl23 — -6. Different values of the mixing angles can shift the normalisation of the 
axes. 

An interesting question is whether we can distinguish between the two or three singlet models by measuring rare 
decay rates? Not unambiguously. If the measured branching ratios disagree with the predictions of the two singlet 
see-saw model, it can be excluded. However, if the branching ratios agree with the two singlet predictions, we cannot 
exclude the three singlet model, which can fit every constellation of branching ratios that is possible in the two singlet 
model. 

Consider now the case where the smallest neutrino mass is small, but non-zero. The Majorana mass eigenvalue M3 
is no longer infinite, but should remain larger than the other two mass eigenvalues Mi^2- However, it is not transparent 
to analytically determine the values of the smallest mass, as a function of other parameters, for which this condition 
is satisfied. A more stringent condition that ensures the validity of our reconstruction is that, in equation (18), for 
each component of M"-*-, the efi^ect of mi is small compared to the effects of the other two mass eigenvalues (normal 
hierarchy) : 

[T>,~^WL\iimi[WlT>z\j « [T>,-^WL]ikmk[WlT>z\j where fc = 2, 3 and /, J = 1, 2, 3 (36) 

Those conditions imply that the reconstruction done in this section is valid if mi << m^, except when Wl is almost 
a diagonal or permutation matrix. In this special case, the contributions of m2,3 are multiplied by a small mixing 
angle or by the cosine of a almost 7r/2 angle, whereas the contribution of mi is multiplied by almost one, so the 
reeconstruction is valid for 

— « GwM or « (Owu - (37) 
mfe V 2 / 

4.2 One zero Yukawa eigenvalue 

In this section, we briefly discuss the case where the lightest neutrino mass vanishes, m^min — > 0, because one Yukawa 
eigenvalue, which we take to be yi, goes to zero. 
Equation (17) can be rewritten: 

WLT>,nWl = T>,VRT>m'^ViT>,vl (38) 



where Dj^ has only two non zero entries on the diagonal and therefore the right hand side of this equation has only a 
2x2 sub matrix with non-zero entries. Furthermore, we know that the doublet neutrino mass eigenstates in do 
not couple to the vanishing first Yukawa eigenvalue jji. Therefore the two angles 6wi2 and 9wi3 are identically zero, 
and a 2 X 2 singlet inverse mass matrix can be reconstructed using eqn (18). However, on the right hand side of 
eqn (38) still contains the three singlet masses. Three equations do not determine three mass eigenvalues and three 
mixing angles, so the whole high energy Lagrangian cannot be reconstructed. 

This limit reduces to the case of two right handed neutrinos. The right hand side of equation (38) can be rewritten 
in terms of two "effective" mass eigenvalues and one "effective" mixing angle. Of course the effective values are 
functions of the three initial mass eigenvalues and the three mixing angles, but wc can only determine the effective 
values. So we cannot observe any difference between the two models unless we can measure the operator [d]. With this 
operator we can distinguish between the two models, but we cannot perform a complete reconstruction, as discussed 
in the following section. 



5 More singlets and the inverse seesaw 

In this section, we consider the prospects for testing, in principle, the seesaw model wc reconstructed. We suppose an 
ideal world where the matrices P, [m] and [d] could be exactly determined. Then from P and [m], one can compute 
d, and compare its elements to their measured values. 

As a first example, consider a model with three singlets, and mi = jji = 0. To be concrete, we take a normal 
hierarchy for the light neutrino masses; a similar reasoning applies for the inverse hierarchy, but in a different two- 
dimensional subspace of doublet flavour space. 

In the singlet basis where Y^^V is diagonal (with eigenvalues {0,2/2,2/3}), we can reconstruct a 2 x 2 submatrix, 
which we call M~^, of the "true" 3x3 matrix 





















[ M3-/ 


Ms2 


^3-3^ J 



(39) 



We imagine this matrix is reconstructed from P and [m], as was discussed in section 4: one obtains the Yukawa 

cigcnvalucis (two non-zero) and the 3x3 rotation matrix Vl from P. This gives the 2x2 matrix Wl, which can be 
substituted into eqn (18) whose left-hand-side is M^^. 

The matrix M~^M*~^, together with Wl and Dj^, can combined to give a ■■prediction" for [d], expressed in the 
neutrino mass basis: 



[d] 



pred 



WjT)^M-^M*-'^T>^WL 



which differs from the "true" expression for [d] = WlTy^M-^M*-'^Ty^WL, 

\ M- 



[d] 
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(40) 



(41) 



which would correspond to the measured value, in our ideal world. It is clear that the reconstructed and true matrices 
differ, and the determinant of the difference is zero. We could repeat this exercise with four singlets, and find that 
the determinant of the difference matrix is non-zero. So in principle, in this best of ideal worlds, it is possible to find 
traces of additional singlets, even when their (Yukawa) interactions with the Standard Model are irrelevant to LFV 
and neutrino masses. 

Notice that since yi = 0, we not can get access to [M~^]ii, and therefore can not determine the mass of the singlet 
TVf. So a complete reconstruction is not possible. This is as expected, because in the limit where [M~-^]ij — >■ O,for 
i 1, the singlet iVf is decoupled from the Standard Model. 

Finally, consider an "inverse seesaw" model, with two singlets Nj of lepton number -1, and two singlets Sj of 
lepton number 1, who have no Yukawa couplings but a small majorana mass matrix /x/j: 



vltHlef, + fH*[Y''*UNi + MjSiNi + ^S-iSj + h.c. 



(42) 



This is a special case of a four-singlet model with only two non-zero eigenvalues in V, and we just argued that a generic 
such four-singlet model could not be reconstructed even in the best of all worlds. Nonetheless, the supersymmetric 
inverse seesaw is reconstructable. From P and [d], which are defined as in the ordinary seesaw in eqns (15) and (8), 
one can obtain V and M. Then from the effective light neutrino mass matrix 



[m] = Y'^M^VM^^Y'' 



T 2 



(43) 



one can obtain [fi] . This could suggest that the supersymmetric inverse seesaw could fit any observations of neutrino 
mass, mixing and non-unitarity, as well as ia ^/37 [8, 30]. An independent determination of the singlet masses(at a 
collider), could test this model. 



Q 



6 Summary and Discussion 



We considered simple toy supersymmetric seesaws (the lightest neutrino was taken massless) with 2,3 and 4 heavy 
singlets, and usually real parameters. We were interested in the properties of a model that allow it to be reconstructed 
in principle, so wc studied the prospects of reconstructing the high-scale parameters of these models from observations 
at or below a TeV. Wc imagine that such observations could give us access to three matrices in doublet lepton flavour 

space (see eqns (15), (5), and (8)) : 

P = vIdIVl , [m] = VM-^V'^vl , [d] = Y''M-iM-i*Y''^' 

It is a straightforward exercise in linear algebra to show that, if V is an invcrtiblc 3x3 matrix, then from any pair 
of the above three matrices, a 3 x 3 matrix M can be obtained. This gives the masses and mixing angles of three 
singlets {Nj}. Therefore, measuring the third of the above three matrices could test the model. 

This programme to test the seesaw suffers from a serious practical problem, which is that P, m and d cannot be 
measured exactly, if at all. We attempted to circumvent this problem by neglecting phases, and by using "reasonable" 
theoretical expectations (such as that the Y" eigenvalues are hierarchical) for unknown parameters. We also take the 
lightest neutrino mass to be zero. 

Then we studied the question of "principle" , whether it is possible to determine the number of singlets participating 
in the seesaw, and reconstruct their properties. 

• Since the lightest neutrino was taken massless, only two singlets are requireed to reproduce the light neutrino 

mass matrix. As shown by Ibarra [5], the seesaw model with two singlets, and two eigenvalues for Y'', is testable. 
In particular, as discussed in section 3 it predicts a discrete set of possible rates for t — >■ £7, given the rate for 
(U — )• ej. If the BR{iJ, — >• 67) is taken at the current bound of 2.4 x 10~^^, then for a normal hierarchy of light 
neutrino masses, one of the BR{t Ij) ^ 10~^ should be visible at a supcr-B factory (for susy parameters as 
described in the Appendix, giving rate expectations as in eqn (56)). In the case of an inverse hierarchy, a similar 
situation could arise, or — unfortunately — both r — )• ^7 decays could be out of the reach of super-B factories. 

To reconstruct the singlet sector of this model, we assume the largest Y^ eigenvalue to be 1, and fix the smaller 
by setting the lighter singlet mass ~ 10^° GeV, for thermal leptogenesis. 

• The three singlet model can have a massless light neutrino in two limits: either one of the singlets is very heavy 
(Mg"^ — > 0), or there are only two non-zero eigenvalues to in Y'^ (so Ui — > 0). 

— The first case, where the 3x3 Y'^ is invertible, is in principle a reconstructable model, allowing to obtain 
the mixing angles and masses of the three singlets. As discussed in section 4 it differs from the two-singlet 
model, in that the three angles of arc; independent. This means that knowing the rate for one Lepton 
Flavour Violating decay, such as — > ej, does not determine the other two. Instead, for hierarchical Y"^ 
eigenvalues, knowing two lepton-flavour-violating branching ratios, allows to predict the third. 
However, as discussed in section 4, this case is not a very realistic limit: if M3 is sufficiently large that it 
contributes negligeably to the light neutrino masses, it is unlikely to be present in loop corrections to the 
sneutrino masses. This means that the Yukawa couplings of A^3 would not contribute to P, which would 
therefore effectively have only two non-zero eigenvalues (see eqn (33)). 

— The second case, where the smallest of the three Yukawa eigenvalues is negligeable, is equivalent to the 
two-singlet case, except that the 2x2 singlet inverse mass matrix that one reconstructs using a 2 x 2 Y" 
matrix in eqn (18), is only a 2 x 2 submatrix of the full inverse mass matrix M~^. As discussed in section 
5, models which all have two non-zero Y"^ eigenvalues, but have two, three or four singlets, are in principle 
distinguishable because they give different predictions for d. However, it is not possible to fully reconstruct 
the singlet sector. This is straightforard to see in the limit where the additional singlets are decoupled, 
because they interact only via the vanishing Yukawa eigenvalue. 

• Four singlet models are in general not reconstructible — that is, the singlet masses cannot be obtained. However, 
if the singlet mass matrix is of the "inverse seesaw" type, it is reconstructible, in principle, from the sneutrino 
mass matrix (P) the light neutrino mass matrix, and the operator d given in eqn (8). 

Almost any observed rates for £a — >■ i/3j can be compatible with the three-or-more singlet supersymmetric seesaw 
[3, 31]. The precise relation between the ^pl rates, and the seesaw parameters depends on various aspects 

and scales of the soft-supersymmctry-breaking mechanism. To circumvent the (not very well delimited) imccrtainties 
resulting from these "soft" assumptions, we estimated the £a — > ^/37 rates using a basic Minimal-Flavour- Violation- 
like parametrisation (eqn (10)), which could give a "lowest order" decription of the supersymmetric seesaw — and of 
several other models. 

The coefficient of the electromagnetic dipole operator (see eqn (9)) is a matrix with indices in doublet and 
singlet charged lepton flavour spaces. We assumed that its flavour-changing elements were proportional to Y^^P, 
where P is a new "spurion" , or matrix in doublet lepton flavour space. A first approximation to the supersymmetric 
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seesaw is to take P = Y'' Y'^''' . We further assumed that the eigenvalues of P were hierarchical, such that the largest 
contributions to flavour changing matrix elements were always mediated by the largest eigenvalue. Finally, we assumed 

that the upper bound on /i — > 57 imposes an approximate zero in P. These assumptions, which could roughly describe 
the three-singlet supersymmetric seesaw, predict a relation between the r — > /i7 and r — > 67 rates (see eqn (34) and 
figure 1). 

The context of this study, was to explore the prospects of "reconstructing" the neutrino mass generation mecha- 
nism from coefflcients of the effective Lagrangian. We focussed on supersymmetric seesaw models, because seesaws 
have many features which facilitate reconstruction. They contain New Physics at tree level, with well-separated mass 
scales. This means that the heavy propagators can be approximated as \/m? — which would be less straightforward, 
for instance, if the New Physics appeared in loops. In addition, all the new interactions and masses of the seesaw are 
flavoured and bilinear. This means that coeflicients in the effective Lagrangian are constructed by matrix multiplica- 
tion, so by combining "suitable" coefficients, one can solve for all the new matrices. This would be less simple, for 
instance, with R-parity violating trilinear couplings. Nonetheless, these nice features of the seesaw are not sufficient to 
ensure reconstructability, because we obviously cannot reconstruct New Physics in the limit where it decouples from 
us. It would be interesting to know what features are neccessary for a reconstructable model. 
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Appendix 

Our model is supersymmetric, so soft supersymmetry-breaking terms for the "sparticles" must be added to the La- 
grangian obtained from the superpotential of eqn (1). These renormalisable terms are interesting for seesaw re- 
construction, because loop processes can imprint traces of the high-scale theory upon them. A "Minimal Flavour 
Violation[20]-like" expectation for the mass-squared matrix of the slcpton doublets {La] could be of the form 



[m^] = I + C^^Y^Y'^T + C^Y^Y 



(44) 



because, in the high-scale theory described by the superpotential of eqn (1), the two spurions linking the 3 and 3 of 
the doublet lepton space are Y^'Y'^^ and Y'^Y''^ In the charged lepton mass eigenstate basis, Y'^Y'^^ generally gives 
off-diagonal elements of [m^], which would induce Lepton Flavour Violating (LFV) processes [32], such as la — >■ ^/37. 

In a model for supersymmetry-breaking and flavour, the soft supersymmctry breaking parameters can be calculated, 
which gives a prediction for the coefficients C™, and for corrections [33] to eqn (44). An estimate for C™ can be 
obtained in gravity-mediated supersymmetry-breaking scenarios, because the soft breaking parameters are present at 
scales ^ M, where loops involving Higgscs(higgsinos) and singlet (s)ncutrinos give small contributions to the doublet 
slepton masses. These are approximately proportional to Y'^Y'^'^, and are resummed by the renormalisation group 
equations. Assuming that, at some high scale Mx ^ M/, the model contains universal soft terms of the form 

£,o/t = mgLl,ia + ao(y^Z" •iJdiT-^+Y'^a/Z" •if„7V=7) + /i.c. , (45) 

then, in the leading log approximation, the slepton doublet soft masses receive contributions [34, 35] of order 



167r2 



(46) 



where a., (i arc indices in the charged lepton mass basis. If one assumes that all the singlets can be decoupled at the 
same scale Mav, and that toq = ao, then one can approximate ^ 



V^log 



'■X 













a/3 







(47) 



which gives 



where we approximate the log 1. 

Lepton flavour violating radiative decays, I, 



ml 
47r2 



proceed via the operator of eqn (9) 



(48) 



(49) 



^Analytic reconstruction is possible without this assumption as shown in section 3.2.3 of [14]. 
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at a rate given by eqn (12) 

BRQc, ^ l0 7) 



rjig ^ If, 7) 



1927r3 



(50) 



For a given weak-scale SUSY spectrum, the coefficients Xp^p can be calculated [36, 37] as a function of masses and 
mixing angles. 

For a = 13 = n, the dipole operator (9) also contributes to the anomalous magnetic moment of the muon 



5-2 



= = ARe{XLnn}mi^ 



(51) 



whose experimental determination [38] differs from the SM prediction by several standard deviations [39]. If this 
anomaly Ja^ ~ 29 x 10~^" [39] is dominantly induced by the same New Physics diagram as £a (-01 (for instance, a 
chargino-slepton loop — see [40] for intuitive estimates of the various diagrams in the mass insertion approximation), 
then [41, 40] 



5a,, 



2 2 



tan /3 



and it is possible to "normalise" [42] the amplitude for — > £^37 by 6a^ 



Xlc 



m 



afj 



5a,i 
4ml 



m 



rric 



TO*' 4TOf, TO^ 

For = TOq, to^^ from eqn (46), and X^afi of the form (10): 
this gives 



and 



BR{lo, ^ Ip 7) 



3Q!|(5a^P 



167r2TO?, 



47rTO4G| 



[Y^Y-+]„^ 



10" 



a0 



(52) 



(53) 



(54) 



(55) 



(56) 



As is well-known [31, 43], the upper bound BR{ji 67) < 2.4 x lO'^^ [22] restricts [Y''Y'''l']^e ^ 2 x 10^^, which 
implies, with the hierarchical assumption for Y'' of eqn (14), that there is a small entry in the mixing matrix Vl'- 

yl\VL^^^Vl^,\ ^ 2x10-3 . (57) 

From current bounds [23] on r ^ /i7 and r — 67, similar approximations give | [Y'^Y'^^]t-£| ^ 0.5. 

Notice that for tan/3 = 10, saturating the [g — 2)^ anomaly via eqn (52) gives fhsuSY ~ 225 GeV, a scale where 
colliders do not seem to find a supersymmetric menagerie. The estimate (56) scales with various parameters as 



10" 



tan;9y /225GeVV 



10 ) \msusY ) 



[Y'^Y'^^] 



(58) 



so changing a parameter by a factor of 2, can change the rare decay rates by an order of magnitude. The "approximate 
zero" of eqn (11), and the expectations which follow from it, relies on the factor in curly brackets in (58) being larger 
than 10-3. 
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